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We report a numerical investigation of tlie fluctuations of the Lyapunov exponent of a two di- 
mensional non-interacting disordered system. While the ratio of the mean to the variance of the 
Lyapunov exponent is not constant, as it is in one dimension, its variation is consistent with the 
single parameter scaling hypothesis. 



I. INTRODUCTION 

The single parameter scaling (SPS) hypothesis ^ is the 
foundation of our understanding of Anderson localization 
in disordered systems. According to the SPS hypothesis 
Anderson localization phenomena are governed by a sin- 
gle parameter: the ratio of the system size to the local- 
ization length ^. When applied to the zero temperature 
conductance g of disordered mesoscopic system the hy- 
pothesis implies that the probability distribution p{g) of 
the conductance obeys^ 



pig)^Fig;L/0. 



(1) 



Here g is in units of 2e^/h and L is the system size. 
The SPS hypothesis has been applied to other physically 
interesting quantities including the localization lengthi^ 
and Lyapunov exponent spectraS of quasi-lD systems, as 
well as the energy level statistics of disordered systems;- 
The probability distribution of all these quantities should 
have a form similar to Eq.(l) if SPS holds. 

Our understanding of scaling is most complete for one- 
dimensional (ID) systems. There are two properties of 
ID systems that distinguish them from higher dimen- 
sional systems. First, their electronic eigenstates are, 
with very few exceptions, always localized. Second, the 
localization length ^ is comparable to the mean free path 
so that there is no diffusive regime. It has been shown 
that, in ID, the cumulants of lug all scale linearly with 
length L.^ It follows that p{g) is log-normal when L ^ ^iS 
The log-normal distribution is determined by two param- 
eters, the mean and variance of Ing. Consistency with 
the SPS hypothesis requires that they be related. For 
weak disorder a perturbative analysis of the ID Ander- 
son model reveals^ 



= ((ln.g-(ln.g))') =2(-lng). 



(2) 



Angular brackets mean an average over disorder. For the 
ID Anderson model, only weak disorder is relevant since 
for strong disorder ^ is comparable to the lattice spacing. 

Equation (j^J holds for many models. The precise con- 
ditions for its validity isiSiiii^ii^ 



e»4 
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where ig is a length scale that is related to the integrated 
density of states. SPS is violated at the boundaries of 
the original spectrum of the system and for fluctuation 
states arising due to disorder in the initial band gaps. A 
violation®'^'* of SPS at the band center of the Anderson 
model was shown in Ref 15 to arise for similar reasons. 

Single parameter scaling of the conductance distribu- 
tion (Eq.(^) has also been verified numerically in the 
three-dimensional (3D) Anderson model close to the crit- 
ical point of the Anderson transitioniiSiil The region of 
validity of the scaling in 3D, however, is not known. One 
can imagine, however, that an inequality similar to Eq. 
may be applicable in this case as well. 

The situation in two dimensional (2D) systems is cur- 
rently very controversial. According to Refi2 all states in 
2D are localized. At the same time there are a large num- 
ber of experiments, in which an apparent metal-insulator 
transition has been observed. (For a recent review see 
Ref 18). The physical meaning of these observation is 
not yet understood, despite a debate that has already 
lasted a decade. The validity, or otherwise, of SPS in 2D 
is, therefore, an important issue. Even for single particle 
models this issue has not yet been fully resolved. For in- 
stance, careful numerical analyses^ii^ of the 2D Anderson 
model showed excellent agreement with SPS. While other 
studies^"'^^'^^'^'^ suggested the existence of power-law lo- 
calized states and two-parameter scaling. Violations of 
SPS have also been reported in more recent papersMi^ 

The example of ID systems demonstrates that con- 
clusive results regarding scaling properties can only be 
obtained from studying the distribution functions of rel- 
evant quantities. Numerical studies^®'^^ in 2D show that 
In g is normally distributed in the regime of strong local- 
ization. It follows that single parameter scaling must be 
manifest in a relation between the average of In g, and 
its variance, similar to the ID equation (PJ- However, at- 
tempts to verify this relation did not reach definite con- 
clusions because of the small system sizes simulated and 
an approach to Eq. Q that was too naive. 

The main objective of our paper is to perform a careful 
analysis of the statistical properties of a 2D disordered 
system of non-interacting electrons, and verify that they 
are consistent with SPS. The object of our calculations is 
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the finite length Lyapunov exponent (LE) for a 2D An- 
derson model with diagonal disorder. (The definition of 
the LE and the meaning of the qualification finite length 
is given below.) For a 2D L x L system with L :s> ^ the 
mean of the LE is equal to the inverse of the localization 
length ^. 

The distribution of conductance has been given special 
attention in the literature because it is directly accessible 
in experiments. However, it should be understood that 
the conductance unavoidably reflects properties not only 
of the system in question, but also properties of the con- 
tacts used to measure it. The Lyapunov exponent, on the 
other hand, is an intrinsic property of the disordered sys- 
tem, containing information about spatial distribution of 
the wave functions, which are ultimately responsible for 
all other properties including conductance. 

We find that the distribution function of the LE is ap- 
proximately normal both when L ^ ^ and when L 3> 
This contrasts with the conductance which exhibits not 
only very strong fluctuations but also a signiflcant change 
in the form of its distribution between the diffusive and 
localized regimes. 

We approach the question of scaling by clarifying the 
relation between the average and variance of the LE that 
is implied by the SPS hypothesis in 2D, and checking 
whether numerical data are consistent with it. We find 
that the relation between the mean and variance is char- 
acterized by a single parameter, namely the ratio of the 
system size to the localization length. Thus, we provide 
convincing evidence that the SPS hypothesis is valid in 
2D disordered systems of non-interacting electrons. 



II. MODEL AND METHOD 

A. The transfer matrix for the Anderson model 

We simulated the two dimensional Anderson model 
with Hamiltonian 



H 



(4) 



The first summation is over all sites on an Lt x L square 
lattice i.e. a system of width Lt and length L. The second 
summation is over all pairs of nearest neighbors. We 
imposed periodic boundary conditions in the transverse 
direction and used a "box" distribution of width W for 
the site energies Ci 



_ / l/W \e,\<Wl2 
\e^\>W/2 



(5) 



Lyapunov exponents arise when the time independent 
Schrodinger equation is expressed as a product of ran- 
dom transfer matrices»22, We divide the system in the 
longitudinal direction into L layers. We form vectors 
'^n of length Lt from the wavefunction amplitudes on 



each layer. For an arbitrary energy E we derive from the 
Schrodinger equation the transfer matrix equation 



= Af„ 



(6) 



The 2Lt X 2Lt transfer matrix Mn relates the wave func- 
tion amplitudes on layer n and n — 1 to those on layers 
n and n + 1. For Eq. 10}, "ifn and M„ are real vectors 
and matrices, respectively, and the transfer matrices are 
identically and independently distributed random matri- 
ces. 



B. Definition of Lyapunov exponents 

We start with a 2Lt x 2Lt orthogonal matrix Qq. We 
perform L transfer matrix multiplications and factor the 
result into a product of an orthgonal matrix Q, a diagonal 
matrix D with positive elements, and an upper triangular 
matrix R with unit diagonal elements 

Ml--- MiQo = QDR. (7) 
We define 2Lt finite length LEs j'^j^^ - - - "f^^"^ by 



(") 



In Ai, 



(8) 



Here D„ is the nth diagonal element of D. The finite 
length LEs are random variables that fluctuate as we 
sample the random potential. For fixed Lt, when L oo 
the LEs always tend to the same limiting values 



r («) 
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(9) 



for (nearly) all samplings of the distribution of transfer 
matrices and (nearly) any choice of Qo^-^ 



The Lf^ LE is the most physically significant: 7 



(Lt) 

is the inverse of the localization length A of an electron 
on an infinite quasi-lD system of width Lt described by 



(10) 



Therefore, in what follows we focus on "/\^*' , dropping 
the superscript and referring to it as the LE 



(Lt) _ 

Tl = 7L- 



(11) 



In numerical calculations, if only the first m LEs are 
required, it is sufficient to make Qq a 2Lt x m real matrix 
with orthonormal columns. Depending on the value of m, 
this can save a considerable amount of computer time. 
The values for the first m LEs obtained in any particular 
calculation are independent of whether or not LEs with 
higher indices are also calculated. 
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To avoid numerical difficulties with the transfer matrix 
multiplication Eq. |(7J), we performed additional Gram- 
Schmidt orthogonalization after every 8 transfer matrix 
multiplications . ^ 

The definition, given in Eq. ©, of the LE for finite 
length that we have adopted here is not the only rea- 
sonable one. We compare our choice with an obvious 
alternative in Appendix ^ 



C. Special considerations for systems of finite 
length 

In this paper we are concerned with the the distribu- 
tion of 7l for finite length L rather than with its asymp- 
totic value as L — + oo. Therefore, we have to deal prop- 
erly with effects related to the finite value of L, effects 
that were routinely considered unimportant in previous 
studies. 

In the asymptotic limit L ^ oo, the value of 7^ de- 
pends only on the distribution of the transfer matrix 
p{Mn), and is independent of the choice of the initial 
matrix Qq. For finite L, however, the distribution of 7^ 
depends on Qq and L, in addition to p(M„). The de- 
pendence on Qq would, if not dealt with, introduce an 
arbitrary element to our analysis that is undesirable. 

To remove the dependence of the distribution of jl 
on Qq we used the following observation to our advan- 
tage. For (almost) any Qq, the distribution of the matrix 
Q approaches an L independent stationary distribution 
Ps (Q) as L increases. The form of Ps (Q) depends only 
on p (Mn). By sampling Qq from ps (Q) we obtain a dis- 
tribution for 7i that depends only on L and p (M„). 

To generate matrices with the required stationary dis- 
tribution Ps (Q), we took an arbitrary set of orthonor- 
mal vectors, performed Nr transfer matrix multiplica- 
tions and factored the result according to Eq. Q. To 
determine how large Nr should be to get a good ap- 
proximation to Ps (Q), we checked whether or not the 
Kolmogorov-Smirnov test could distinguish between the 
distributions of for i = 1 obtained with different Nr. 
The test showed that once Nr > 100, the distribution 
of the LE becomes independent of N^. Below we set 

Nr = 1000. 



III. RESULTS 

Since our interest in this paper is in the distribution of 
the LE in 2D systems we set the width Lt and length L 
of the system equal 

Lt = L (12) 

i.e. in the remainder of the paper we consider only 2D 
L X L systems. 



A. Distribution of the LE 

We simulated systems with Fermi energy E = 1, dis- 
order 5 < W < 14 and a range of systems sizes between 
L — 16 and L = 512. The distribution of the LE for two 
particular cases are shown in Fig. ^and Fig. [3 These are 
representative of the parameter range we studied. Figure 
n corresponds to the situation L ^ ^, while Fig. |21 cor- 
responds to the situation L ^ ^. In the figure captions 
we give the values of the mean, variance and skewness 
for the numerical data, as well as the number of samples 
simulated. 

The skewness is a measure of the symmetry of the dis- 
tribution. Distributions that are symmetrical about their 
mean, such as the normal distribution, have a skewness 
equal to zero. According to Ref. [s^, a distribution whose 
skewness has absolute value greater than unity is con- 
sidered highly skew. A distribution whose skewness has 
absolute value less than one half is considered fairly sym- 
metrical. For data sampled from a normal distribution, 
the skewness is expected to be distributed around zero 
with a standard deviation of y^l5/Ns where Ng is the 
number of samples^ 

For the data in Fig. Q the difference of the skewness 
from zero is not statistically significant. This is consistent 
with the LE having a normal distribution. 

For the data in Fig. El the difference of the skewness 
from zero is statistically significant. What is the physical 
significance of this deviation? Normally we would expect 
the scaling hypothesis to apply only when the localization 
length ^ is much longer than microscopic length scales 
such as the mean free path, lattice constant etc. Here 
these are approximately unity, so this condition corre- 
sponds to f ^ 1. This condition is satisfied for the data 
in Fig. n where ^ ~ 180 (see Table Hj), but not for the 
data in Fig. [21 where ^ ~ 2.5. Therefore, we think that 
the deviation from the normal distribution seen in Fig. 
121 is not significant in the context of our study. 

In our opinion, the normal distribution is a reasonable 
approximation to the observed distribution for the range 
of i/^ in our simulations. In what follows, we concen- 
trate our attention on the mean and variance of the LE 
and their scaling. 

It is also important to bear in mind when looking at 
Figs, n and that the scaling hypothesis is expected to 
apply to the bulk of the distribution not its tails, i.e. to 
typical states not necessarily to very rare states. Hence, 
we use a linear scale for the probability density axis and 
not a logarithmic scale, which would unduly emphasize 
the tails of the distribution. 



B. The scaling of the mean LE 

According to the SPS hypothesis the scaling of both 
the mean and variance of the LE should be governed by 
the same length scale, the localization length ^. A quan- 
titative test of scaling involves checking the consistency 
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FIG. 1: The distribution of the LE for a 2D L x L system 
with E — 1, W = 5 and L = 32. The hne is a normal distri- 
bution with mean and variance equal to that of the numerical 
data. The numerical data have mean 0.032, standard devia- 
tion 0.032 and skewness -0.0031. The number of samples is 
65,523. 
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FIG. 2: The distribution of the LE for a 2D L x L system with 
E = 1, W = 14: and L = 128. The line is a normal distribution 
with mean and variance equal to that of the numerical data. 
The numerical data have mean 0.40, standard deviation 0.035 
and skewness -0.33. The number of samples is 40,000. 



of the disorder dependence of ^ obtained independently 
from the scaling of the mean and variance of the LE. In 
this section we deal with the scaling of the mean LE. 

For the scaling analysis of the mean LE we esti- 
mated (7l) to a precision of 0.25% for system sizes 
L = 16,32,64,128,256 and 512. For W < 11 the max- 
imum system size was reduced to 256. We determined 
£, = £, {W) by fitting the variation of the mean LE with 
L and W to the SPS law 



(7l) L^F 



(13) 



When L ^ ^ we suppose that the mean of the LE will 
tend to the inverse of the 2D localization length ^ i.e. 



lim (7/ 



This is equivalent to 



F{x)^x a;>l. 



(14) 



(15) 



For numerical reasons we expressed the scaling function 
in the form 



logio {il)L^ f [ log 



(16) 



and used a spline to interpolate the function /. The val- 
ues of /, at the values of L/^ in Tabled were fitting pa- 
rameters. To ensure the spline interpolation reproduces 
Eq. 113, we fixed the value of / at L/C = 1000. The 
corresponding value of F is given in parenthesis in Table 
^ The remaining fitting parameters were the localization 
lengths for each disorder. Finally, we used the shape pre- 
serving Akima spline to avoid unphysical oscillations of 
/. We summarize the results in Tables and and in 
Figure ini We can see from this figure that the data for 
different values of disorder W and system size L fall on 
a common scaling curve when expressed as a function of 
L/^. Moreover, for large L we observe the expected linear 
dependence, with slope equal to the inverse localization 
length. 



C. Scaling of the fluctuations of the LE 

Taking into account that the dimension of 7^ is 1/L, 
and that of its variance cr^ is 1/L'^, we can define a di- 
mensionless quantity S by 



E = 



(7l)' 



(17) 



According to the SPS hypothesis the localization length 
is the only relevant length in the system, so S should 
obey the SPS law 



(18) 
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TABLE L The 2D localization length and the scaling function 
F determined from the scaling of the mean LE. The errors 
quoted are 95% confidence intervals and are estimated using 
the Monte Carlo method. 



w 






F 


5 


178 ±2 


0.1 


0.86 ± .005 


5.5 


85 ±1 


0.5 


1.51 ± .01 


6 


48 ±.5 


1 


2.11 ± .02 


6.5 


30.2 ± .3 


2 


3.17 ± .02 


7 


20.8 ± .2 


5 


6.16 ± .05 


8 


11.7± .1 


10 


11.0 ± .1 


9 


7.69 ± .08 


20 


20.6 ± .2 


10 


5.54 ± .06 


1000 


(1000) 


11 


4.26 ± .04 






12 


3.44 ± .03 






13 


2.88 ± .03 






14 


2.47 ± .02 







TABLE II: Details of the finite size scaling fits: the number 
of data Nd, the number of parameters Np, the value of for 
the best fit and the goodness of fit probability Q. 



Statistic and data range 


Np Na Q 


16 < L < 512, 5 < W < 14 


19 63 46.0 0.4 


S 

64 < L < 256, 5 <W <12 


15 30 20.3 0.2 



100- 



10- 



1 : 




1 10 100 



FIG. 3: The finite size scaling fit (line) to the data (circles) 
for the mean of the Lyapunov exponent. The precision of the 
numerical data is 0.25%. 



In ID, the linear scaling of the cumulants of Ing, and 
the relation between the LE and g described in the Ap- 
pendix, allow us to deduce from Eq. (jSJ the much more 
prescriptive statement 

E = l. (19) 

However, for a 2D L x L system the cumulants of In g do 
not scale linearly with L, except perhaps in the regime 
where _L 3> ^; a regime which it is more difficult to reach 
in 2D than in ID. Therefore, we should not expect that S 
be unity or even constant in our calculations. Confirma- 
tion of the SPS hypothesis in 2D consists not in demon- 
strating that calculated values vary in accord with Eq. 
(|19|l but rather in trying to establish Eq. H18|) . 

We have plotted the variation of E with system size 
in Figure 01 where different lines correspond to different 
values of disorder, W. These data were analyzed in an 
analogous way to the mean LE. We expressed the SPS 
law Eq. (fl^ in the form 

E = /s (^logio 1^ , (20) 

and used an Akima spline interpolation of the function 
/s. The values of /s, at the values of L/^ hsted in Table 
mil and the localization lengths for each disorder were 
fitting parameters. To obtain a reasonable goodness of 
fit (> 0.1) we had to restrict the range of data considered 
to 5 < < 12 and L > 64. (There seem to be a 
more pronounced finite size correction in the data for 
the variance than in the data for the mean LE. Also, 
the breakdown of scaling when ^ is comparable to the 
lattice spacing (^ ~ 1) seems to be evident sooner in the 
variance of the LE than in the mean LE.) 

When fitting data for the mean LE, we were able to 
determine the absolute value of ^ with the aid of H15|l . 
Unfortunately, no similar relation is available for /s and 
so we cannot fix the absolute scale of ^ by fitting data 
for E alone. Indeed, looking at Fig. Owe can see that, 
if we translate both the fit and the data by the same 
amount parallel to the abscissa, we obtain an equally 
good fit. To avoid this ambiguity, we set the value of 
the localization length for = 12 to that found for 
the mean LE. We show the results in Table lllll Apart 
from an over estimate of ^ for W = 5 and W — 5.5, 
the results are consistent with those for the mean LE. 
In addition, in Figure |5l we have plotted E versus L/^, 
where ^ is estimated from this fit. We see that all the 
different curves of Fig^ collapse on to a single curve, 
confirming the correctness of Eq. (|18|l . We conclude 
that the fluctuations of the LE are consistent with the 
SPS hypothesis. 

Looking at Fig. O it is plausible that the function E 
will tend to a finite asymptotic value as i/^ ^ oo. If this 
does occur, the fluctuations of the LE in the 2D asymp- 
totic limit (L/^ ^ oo with Lt = L) decay as 1/VX. This 
is similar to the behavior in the quasi-lD limit (L/^ — > oo 
with Lt flxed) where the fluctuations in the LE also decay 
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FIG. 4: E versus system size. For each point 160, 000 samples 
were simulated, corresponding to a precision of approximately 
0.4%. The lines, which connect points corresponding to a 
common value of the disorder W , are a guide to the eye only. 



TABLE 111: The 2D localization length and the scaling func- 
tion F determined from the scaling of E. 



w 






Fe 


5 


226 ± 38 


0.1 


.98 ± .02 


5.5 


96 ±8 


1 


.948 ± .005 


6 


50 ±2 


5 


.804 ± .004 


6.5 


30.4 ± 1 


10 


.683 ± .004 


7 


20.9 ± .6 


20 


.551 ± .003 


8 


11.9 ± .3 


50 


.403 ± .002 


9 


7.8 ± .2 






10 


5.6 ±.1 






11 


4.3 ± .1 






12 


(3.44) 







as 1/ ^/L. The only difference is that E is always unity 
in quasi-lD, while the asymptotic value of E is less than 
unity in 2D. 



IV. CONCLUSION 

We have investigated numerically the scaling of the 
fluctuations of the LE in the 2D Anderson model with 
diagonal disorder. We found that the distribution of the 
LE is approximately normal both when L <C ^ and L ^ 
5. We showed that the parameters of the distribution, 
the mean LE and its variance behave in accordance with 



1.0 -■ 



0.8 - 



0.6 - 
0.4 - 



0.2 - 



0.0 — • — • — I 

1 10 100 

FIG. 5: A visual check that the numerical data for E defined 
in Eq. ^ (points) satisfy the SPS law Eq. lHHJ (line). 



the single parameter scaling hypothesis for the energy 
value considered in our calculations. This value, E = 1, 
was chosen to lie far from the boundaries of the initial 
spectrum, E = ±4, and from its center E = in order to 
avoid anomalies related to the band edge behavior, which 
were found in ID systems. We expect that the behavior 
for any other value of the energy will be similar as long as 
it is not close to an anomalous region. We found that the 
manifestation of SPS in numerical studies of 2D systems 
is different from that of ID systems. Instead of the simple 
relation between the mean and variance of the LE given 
by Eq. p9(l and valid for ID, in 2D one has to analyze 
the compliance of numerical data with the SPS relation 
Eq. (CHI). 

The fact that we verified SPS both when L ^ and 
i 3> ^ is smiificant because it contradicts the conclu- 
sions of Ref^3 and Ref f25l where a behavior inconsistent 
with SPS was found. While a complete elucidation of the 
sources of this disagreement is beyond the scope of this 
paper, we can suggest some possibilities that might be 
worth pursuing in future work. First, it is possible that 
the logarithmic increase of the localization length seen 
in Rcf^i! might be reconcilable with the SPS hypothesis, 
in much the same way as we have reconciled the system 
size and disorder dependence of the ratio E with SPS 
here. Second, the authors of RefS and Ref'25' analyzed 
the spatial properties of eigenfunctions. In ID the rela- 
tionship between lengths that characterize transport and 
wavefunctions is well establishedi^ In 2D there maybe 
aspects of this relationship that have not yet been prop- 
erly understood. Third, in R,efl23 and R,efl25l wavefunc- 
tions corresponding to i? = were studied. In ID this is 
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a special spectral point, at which SPS is violated j"'^ i It 
seems reasonable to suggest that E — is also a special 
point for 2D where SPS should not be expected. 

APPENDIX A: ALTERNATIVE DEFINITION OF 
THE LE FOR FINITE LENGTH. 

The definition of the finite length LE we have used in 
this work is not the only reasonable one. In this appendix 
we will describe an alternative and compare with the def- 
inition described in the main text of this paper. Given a 
transfer matrix M 

L 

M = [] M„, (Al) 

we can define a matrix by 

n = liiM^M. (A2) 

The eigenvalue spectrum of Q is composed of pairs of 
opposite sign — : n = I ■ ■ ■ Lt} . From these 

eigenvalues we could define the LE in an alternative way 
as 



In the limit that L — > oo at fixed Lt, the random vari- 
ables defined by Eq. ljA3|l always tend to the same lim- 
iting values for (nearly) all samplings of the distribution 
of transfer matrices. These values are the same as those 
obtained with Eq. (O in the same limit. For finite L the 
values of Eq. © and Eq. (|A3p are different. We sum- 
maries the main characteristics of each definition below. 
For the definition Eq. ^ in the main text: 

PI The LE are not the eigenvalues of a matrix. The 
indices of the LE refer to the order in which they 
are obtained from the Gram-Schmidt procedure. In 
general, this is not in a strictly decreasing order. 

P2 Though the sum of the all LEs is always zero, for 
finite L and for a single sample, the LE do not occur 
in pairs of opposite sign. This symmetry is restored 
after taking the limit L — > cxd for a single sample, or 
after averaging over an ensemble of samples. For a 
single sample we have found that the symmetry also 
appears when Qq is sampled from the stationary 
distribution's (Qo) described in III Gl 

P3 For fixed Lt and Qo sampled from ps (Qo), the mean 
of the LEs are independent of L. (Note that in the 
main text we consider scaling with Lt = L, so this 
property is not applicable there.) 



P4 The LE have a simple geometrical interpretation in 
terms of the exponential rate of increase of lengths, 
areas, volumes etc. 
For the definition (|A3|) : 

PI The LEs are related to the eigenvalues of a matrix 
and hence there is no prescribed ordering for them. 
It is conventional to put the LEs in decreasing order 
and the index in the definition (|A3I) usually refers 
to this order. 

P2 For all L, the LE occur in pairs of opposite sign. 
This exact symmetry is exhibited not just after av- 
eraging over an ensemble of samples but also by a 
single sample. 

P3 In this definition there is no analogue of Qo smd 
hence no analogue of property P3 for Eq. ©. 

P4 There is no simple geometric interpretation except 
in the asymptotic limit. 

SPS can be investigated using either definition. The 
quantity defined by Eq. © has the advantage that its 
distribution is normal, while at the same time retaining 
a straightforward relationship to the decay of the wave- 
function in the disordered system. 

APPENDIX B: RELATION OF LE TO 
CONDUCTANCE IN ID. 

For a strictly ID system whose length is much longer 
than the localization length the transmission coefficient t 
for the transmission of electrons through the disordered 
sample decays as 

-ln|i| =lnDi +0(L") (L > ^ • (Bl) 

The 0{L^) term is a fluctuating term that depends on the 
nature of the leads attached to the sample when defining 
the scattering problem. Using the Landaucr formula to 
relate the transmission and the conductance we have 

^\ng = 2jLL + 0{L") (L > ^ . (B2) 

From this we deduce that equation Eq. ^ is equivalent 
to equation Eq. H19I) when L ^ ^. 
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